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Abstract 

The fluctuation-dissipation theorem (FDT) is very general and applies to a broad variety 
of different physical phenomena in condensed matter physics. With the help of the FDT 
and following the famous work of Caldeira and Leggett, we show that, whenever linear 
response theory applies, any generic bosonic or fermionic system at flnite temperature 
T can be mapped onto a flctitious system of free harmonic oscillators. To the best of 
our knowledge, this is the flrst time that such a mapping is explicitly worked out. This 
flnding provides further theoretical support to the phenomenological harmonic oscillator 
models commonly used in condensed matter. Moreover, our result helps in clarifying an 
interpretation issue related to the presence and physical origin of the Bose-Einstein factor 
in the FDT. 

1 Introduction 

The idea of modeling physical systems as a collection of harmonic oscillators has a long 
history and dates back to even before the birth of quantum mechanics. One of the best 
known example is Planck's work on black body radiation at the edge of the classical 
(beginning of quantum) era. More recently, this model has been of great importance in 
connection with the study of dissipation in quantum mechanics [H [21 E]- 
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In their famous paper devoted to the study of tunnehng in dissipative systems [1], 
Caldeira and Legget observed that "any physical system which is weakly perturbed around 
its equilibrium state can be adequately represented (at T=0 at least) by regarding that 
system as equivalent to a set of simple harmonic oscillators". They supplemented this 
statement with an explicit computation, where they showed that given a quantum system 
at T = 0, in the lowest order approximation, its dynamics can be reproduced with the 
help of a properly constructed system of harmonic oscillators. They also stressed that the 
study of the T 7^ case (not considered in their paper) needed separate discussion. 

Caldeira and Legget then suggested the form of the total lagrangian of a physical 
system in interaction with a certain "environment" (the "heat bath") and worked out 
the consequences of this assumption in connection with the tunneling problem [H 
Prompted by this pioneering work, harmonic oscillator models are nowadays extensively 
used and there is little doubt that, from a phenomenological point of view, they reproduce 
quite well the physics of the systems under investigation. As a specific example, we 
can consider an electrical circuit, where the resistance is modeled with a collection of 
harmonically oscillating electrical dipoles. 

From a theoretical point of view, however, it would be more satisfactory if we could 
prove that, for any generic system at finite temperature T 7^ 0, it is possible to find 
an equivalent system of harmonic oscillators such that the statistical (thermodynamical) 
properties of the real physical system are properly reproduced by the system of oscillators. 
This would be an extension of the Caldeira Legget result and would provide further 
theoretical support to the commonly used phenomenological models. 

The main scope of this work is to present a new and very general result which pro- 
vides the above mentioned extension of the Caldeira-Legget one. By working within the 
framework of the FDT[6], we show that, whenever linear response theory applies, any 
generic bosonic and/or fermionic system at finite temperature can be mapped onto the 
Fock space of a fictitious system of free harmonic oscillators at the same temperature. 

As a byproduct of our analysis, we shall see that our finding should help in clarifying an 
interpretation issue concerning the Bose-Einstein (BE) distribution factor which appears 
in the FDT. Actually, an often raised question concerns the physical meaning and/or 
origin of the BE factor which appears in the relation between the power spectrum of 
the fluctuating quantity and the corresponding generalized susceptibility. Sometimes this 
term is interpreted as due to an harmonic oscillator composition of the physical system 
under investigation. Such an interpretation, however, is not supported by the derivation 
of the theorem itself (see for instance [3 El E]). Moreover, the FDT applies to any generic 
bosonic or fermionic system (irrespectively of its statistics). 

Far from being an academic question, the resolution of this interpretation issue is 
of very practical importance in many different contexts [TT| [T^ [T5] . From a real 
understanding of the origin of this term often depends the correct physical interpretation 
of theoretical and experimental results [H [9l [TU [151 IIHI HZ]- As we shall see, our results 
suggest that this term does not originate from underlying physical oscillator degrees of 
freedom of the system but is rather a general property related to the approximation (linear 
response) involved in the derivation of the FDT. 
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According to this theorem, whenever hnear response theory is apphcable, given a 
generic system which interacts with an external field f{t) through the interaction term 
V = —f{t) A, where A is an observable of the system, the mean square of the Fourier trans- 
form A{u) of A{t) is related to the imaginary part x'^i^) '^^ corresponding (Fourier 
transformed) generalized susceptibility by the relation [B] : 

{A'{.)) = fi,:(.)coth (■^'^ = 2fi,:(.) (i + . (1) 

where /3 = 1/kT, with T the temperature of the system and k the Boltzmann constant. 

For instance, in the case of a resistively shunted Josephson junction [10], when applied 
to the power spectrum 5*7(0;) of the noise current (fiuctuation) in the resistive shunt 
(dissipation), the theorem takes the form {R is the shunt resistance) 

S,M = ^(- + j5^j. (2) 

The power spectrum Sj{u!) has been measured |1U) and good agreement between the ex- 
perimental results and Eq. ([2]) has been found. 

The above ^ term is sometimes presented |10l [TT| [T8| [T9] as due to zero point energies 
and the experimental results |10) as a measurement of them. In fact, the term in paren- 
thesis in Eq. ([2]) coincides with the mean energy of an harmonic oscillator of frequency u 
in a thermal bath. The same holds true for the general case of Eq. ([T]), where the similar 
term is the mean energy of an harmonic oscillator in hw units, i.e. the BE distribution 
function. In the following we shall see that our result (the mapping) strongly suggests 
that the agreement between the experimental results |10| and Eq. ([2]) cannot be considered 
as a signature of measurement of zero point energies. 

The rest of the paper is organized as follows. In section 2 we briefiy review the 
derivation of the FDT and establish some relations useful for the following. In section 
3 we establish our new result, the mapping, i.e. we show that for any given bosonic or 
fermionic system at finite temperature T 7^ we can always find a fictitious system of 
harmonic oscillators in such a manner that the physical quantities which appear in the 
FDT can be obtained from this equivalent system of oscillators. Section 4 is for our 
comments and conclusions. In particular, in this last section we present our comments on 
the interpretation issue related to the BE term in the FDT. 



2 The fluctuation-dissipation theorem 

Let us begin by briefiy reviewing the derivation of the FDT. Consider a macroscopic 
system with unperturbed Hamiltonian Hq under the infiuence of the perturbation 

V = -mA{t), (3) 
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where A{t) is an observable (a bosonic operator) of the system and f{t) an external gener- 
alized forc^. Let \E„) be the eigenstates (with eigenvalues £"„) and {En\A{t)\En) = 0. 
Within the framework of linear response theory, the quantum-statistical average {A{t)) f 
of the observable A{t) in the presence of V is given by 

{Mt))f= f dt'xAt-t')m (4) 
where Xa(^ ~ generalized susceptibility, 

xAt - 1') = \e{t - = -\GK{t - 1') , (5) 

with (...) = ZnQn{En\...\En), Qn = e"^^"/;^ , Z = En e"^""" . ^^(t - t') being the 
retarded Green's function and A{t) = Q-^^ot/hy^^-iHot/h^ 
Defining the correlators (from now on t =0): 

G>(t) = {A{t) i(0)) and G<(t) = (i(0) A{t)) , (6) 

so that Gait) = —i6{t){Gy{t) — G^{t)), and the corresponding Fourier transforms, G^{u) 
and G^{uj) respectively, it is a matter of few lines to show that: 

G^uj) = --^-^ImGniu) ; G<(a;) = e"^''- G>(a;) . (7) 
Finally, by noting that 

(i'M) = ^(G>(u;) + G<(u;)) (8) 

and that the Fourier transform of XaI'^) Xa('^) ~ ^1^'^) ~'~ ^x'^i^) ~ ^jfiRi^) S^t: 

{A\u)) = hxM^^^ = hxJu;)coth^^j=2hxAu^) (2 + ;^^; 

which is Eq. ([T]), the celebrated FDT. 

As observed by Kubo et al. [7j (and shown in the derivation sketched above), the BE 
factor is simply due to a peculiar combination of Boltzmann factors in Eq. ([9]) and there 
is no reference to physical harmonic oscillators of the system whatsoever. Despite such an 

authoritative remark, some people insist in interpreting the ^| + ^phL_i ^ term as related 
to harmonic oscillator degrees of freedom of the physical system. 

In the case of the measured |10] power spectrum of Eq. (|2]), some authors [T2| [T3] 
interpret this term as due to the electromagnetic field in the resistive shunt and therefore 

"^More generally, we could consider a local observable and a local generalized force, in which case 
we would have V = — J cfi fA{r)f{r,t), and successively define a local susceptibility '''j ^') (see 

Eq. ^ below). As this would add nothing to our argument, we shall restrict ourselves to r-independent 
quantities. The extension to include local operators is immediate. 
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the first term in parenthesis of Eq. ([2]) as originating from zero point energies of this 
electromagnetic field. Such an interpretation, however, is not supported by any physical 
derivation and has been strongly criticized in [HI [HI [151 [IS]- Very recently, starting from 
the results of the present work, we have also carefully investigated this issue [9], providing 
arguments which strongly support previous criticisms (see j^] for details). 

Let us go back now to our analysis. For our purposes, it is useful to show that from 
Eqs. ([5]) and ([8]) we can easily derive the following expressions for Xa{^) ^ind (yl^(cj)) : 
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and 
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Ei — Ej 

h 



+ uj] +6 



h 



Ej — Ei 



h 



(10) 



fill 



where = {Ei\A\Ej) . 

In fact, by inserting in Eq. ([5]) the expressions: 
9{t -t') = - ^^^^ , / = and A{t) = e^^^^m Ae-^^-^l^^^ we get: 



fi j In u + irj 



^Ej){t~t')/n _ -i{E,-E,)(t^t')lh 



:i2) 



Then, making use of lim^ 



As for Eq. (Ilip . from Eq. ([2]) for G^{t), we have: 



V (^) — in6{uj) , Eq. ([TU]) follows immediately. 



G>{t) = J2pmei^'Ae-i^'\E,){E,\A\Ei) = = J^P^e-^^^^-^^'lA,^ ■ 
Working out the similar expression for G<(t), for the correlation function G{t) we get: 
G{t) = he At) + G<(t)) = I J2 P^l A.f (e-^^^^-^'^* + e-i(^-^^)*) , (13) 



SO that the Fourier transform G{u!) is: 
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Ej — Ei 

h 



+ uj] + 6 



Ei — Ej 

h 



(14) 



As = G{0), G{u) is the spectral density {A\u)) of (i^) ( (^2^ = (^A\u)) '1^ ) . 
Then, making use of Eq. f ll4p . we finally get Eq. ([TT]) . 

For our purposes, it is also useful to write Eq. f[TT]) in a different manner. After some 
straightforward manipulations, Eq. ([TT]) can be written as: 



{A\u)) = nJ2iQ^ - Q3)\Aj\' coth (^)[S {CO -uj,.) + 6 {u; + CO,.)] (15) 
j>i ^ ^ 

= 71 coth (^^^^ ^{Qi - gj)\Aij\^ [6 {u - ujji) - 6 {u + Uji)] , (16) 



where we have introduced the notation: ujji = ^ \ By following similar steps, Eq. ( ITOl) 
can also be written as: 

X"i(^) = \ Y^{Q^ - Q'j)\Aj? - + ^3^)\ ■ (17) 

j>i 

Clearly, comparing Eq. f|T6|) with Eq. f|T7|) . we find, as we should, the FDT theorem. 

Now, starting from Eqs. (ITSl) and ( IT71) and taking inspiration from the seminal work 
of Caldeira and Leggett |1] , we shall be able to establish a formal mapping between the 
real system considered so far and a system of fictitious harmonic oscillators. A similar 
mapping, restricted however to the T = case, was considered in where it was also 
noted that the T 7^ case needs separate discussion. The mapping that we are going to 
construct in the present work deals with the T ^ general case. 



3 Mapping Boson and Fermion systems onto harmonic 
oscillators 

To prepare the basis for the construction of this mapping, let us consider first a real 
system Sosc of harmonic oscillators (each of which is labeled below by the double index 
{ji} for reasons that will become clear in the following) whose free Hamiltonian is: 

j>i ^ J* 

where Uji are the proper frequencies of the individual harmonic oscillators and M,j their 
masses. Let \nji) [uji = 0, 1,2, ...) be the occupation number states of the {ji} oscillator 
out of which the Fock space of Sosc is built up. Let us consider also Sosc in interaction 
with an external system through the one-particle operator: 

Vosc=-fit)Aosc, (19) 

with 

Aosc = ^ {oiji qji) . (20) 

j>i 

Obviously, the FDT applied to Sosc gives {Al^^{uj)) = hXosd^) coth (^^) , but this is not 
what matters to us. 

What is important for our purposes is that, differently from any other generic system, 
for Sosc we can exactly compute {A^g^uj)) and Xosd^) from Eqs. (fTOj) and ( ITT]) because 
we can explicitly compute the matrix elements of Aosc- 

In fact, if we apply Eqs. (I13p and f ll4p to Sosc and replace the double index notation 
[uji ; riji ; Mji ; etc.) with the more convenient (for the time being) and self explanatory 
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one index notation (wi, W2, •■■ ; ni, 77-2, ... ; Mi, M2, ... ; etc.), for {Alg^{uj)) we have: 

f?ni f?n2 ■)\{ni,n2, ..\A 

osc |mi,m2, ..;| 

711, 712, .■ mi,m2,.- 

X [5 (w + /iCJi + I2CO2 H ) + 5 (w - /iCJi - /2a;2 )] , (21) 

where h = Uk - rrik, Qn^ = e-/5("fe+V2)na;,/^^^ = J2n^ ^-^inu+imru^k (note also that in 
this one index notation Aosc is written as Aosc = J2k i^k %) )■ Now, as 



{nk\qk\mk) 



we immediately get: 



h 



2MkUJk 



(y/uk + l(nfc + l|mfc) + y/n^{nk - l\mk)) 



(22) 



(AIM) 



ni,n2,.. mi,m2,.. 



X 



E E 

,1,7^ 

E«fc 



2MkUk 



{y/nk + l^mfe,nfc+l + V^'^mfe,nfc-l) H 

X [5 (w + /iwi + I2CU2 H ) + 5 (w - hiui - I2UJ2 -•■■)] • 



(23) 



Let us concentrate our attention to the square in the second hne of Eq. (123|) . Due 
to the presence of the Kronecker deltas, all the crossed terms in this square, i.e. all the 
terms with different values of the index k, vanish. In other words, the square of the sum 
is equal to the sum of the squares: 



E 



ak 



h 



2MkUJk 



{y/Uk + 15mfe,nfe+l + VnkS.. 



'7T!.fc,rife-l 



h^k 



E 



Oik 



h 



2MkUk 



(a/ Uk + 15mfc,nfc+l + \/nkSmk,nk-l) Y\, ^"^h^'^h 1 

hj^k / 



(24) 



k \ ' hf^k 

For the same reason, the same holds true for each value of the index k, i.e.: 



h 



2MkUk 



7Tl J:, 71^.-1 



h 



2MkUk 

Therefore, as h = Uk — rrik, for (Al^^^u)) we get: 



^7nj. ,71^.-1) ^mh,nh • 



(25) 



ni,n2, 



h 

2MkUk 



{2nk + l){6{u-uJk)+S{u + Uk)). (26) 



Finally, as Yln^ ^"■k ~ above expression becomes: 



2MkUk 

k nfc 



E ^M^^ (^^) - a;,) + 5{u + u,)) . (28) 



TT 



Going back to the original double index notation: 

(^LM) = vr E coth V5(c. - c.,,) + 5(c. + c.,^ (29) 

= TT coth f ^ j XI 2M~^ ~ ~ + ^J^)) ■ '^^^^ 

We have just seen that given a real system Sosc of harmonic oscillators and the one 
particle operator Aosc of Eq. (120 p . for such an operator is possible to evaluate explicitly 
{Alg^{uj)). We find that each of the individual harmonic oscillators gives rise to a term 

coth which in turn comes from the term Xln 0nji {'^^ji + 1) of Eq. (I27p . 

Let us now consider Xosc{^)i which (see Eqs. (1101) and (123|) ) is nothing but: 



X [5 (w + /iwi + /2W2 H )- 5 (cj - /iwi - /2t^2 -■■■)] ■ (31) 

Apart from the factor 1/^, Eq. fl3T]) differs from Eq. (|23|) because it contains the difference 
(rather than the sum) of delta functions in the last line. 

If we proceed for XosSm^) have just done for (A^^^(a;)), we immediately note that 

the only difference with the previous computation is due to this minus sign. In fact, its 
presence causes that rather than the combination (2 Uk + 1) of Eq. (126|) . which comes from 
the sum (n^ + 1) + rifc of Eq. (125 p . we get the combination (n^ + 1) — = 1. Therefore, 
for Xosc('^) we do not get the sum Qn^. (2 Uk + l) = coth ('^^) of Eq. ( !27|l . but rather 
En, =1. Then: 



Naturally, comparing Eq. (|30|) with Eq. (1320 we see that for Sosc the FDT holds true, as 
it should. However, what is important for our purposes is to note that for this system we 
have been able to compute separately {Alg^{uj)) and x'osd^) ^^'^ found that the coth (^^) 

8 



factor of the FDT originates from the individual contributions coth y j of each of the 

harmonic oscillators of Sosc- 

We are now in the position to build up our mapping. Let us consider the original 
system S, described by the unperturbed Hamiltonian Hq , in interaction with an external 
field f{t) through the interaction term V = —f{t) A (see Eq. ([2])), and construct a fictitious 
system of harmonic oscillators Sosc, described by the free Hamiltonian Hqsc of Eq- fll8p . in 
interaction with the same external field f{t) through the interaction term Vosc of Eq. ( I19p . 
with Aosc given by Eq. (120|) . where for aji we choose: 

c^ji={ {Qi - Qj)-^ (33) 

and for the proper frequencies Uji of the oscillators: 

Uji = {E,-E,)/h>0, (34) 

with Ei the eigenvalues of the Hamiltonian Hq of the real system. 

Comparing Eq. f lHU]) with Eq. flTB]) and Eq. fl32p with Eq. (1170 . it is immediate to see 
that with the above choices of aji and Uji we have: 

(i2(a;)) = (iLM) and x"M = x'U^) . (35) 

Eqs. fl33|) and are the central results of our analysis. Actually, these are the 
equations which allow to establish our mapping. In fact, with such a choice of the a's 
and the tu's, we are able to map the real system S onto a fictitious system of harmonic 

oscillators Sosc , 

S Sosc, (36) 

in such a manner that x'^i^) {^"^i^)) of the real system are equivalently obtained by 
computing the corresponding quantities of the fictitious one (Eqs. (135|) ). 

This is the desired result. What we have just shown is that any generic boson or 
fermion system at finite temperature T is equivalent to a system of harmonic oscillators 
at the same temperature. From a theoretical point of view, the relevance of such a 
result should be immediately clear. As we have already observed, in fact, harmonic 
oscillator models are quite common in modeling generic physical systems. Now, the 
typical physical situation we have to deal with is that of a system (bosonic or fermionic) 
at finite temperature T ^ 0. In this respect, our mapping fills up the gap mentioned by 
Caldeira and Legget (see Appendix C of [1]) by extending the T = mapping put forward 
by them to the general finite temperature case, thus providing further theoretical support 
to the use of these models. 

In the following section, we would like to add some more comments on the above 
results. In particular, we are going to consider the previously mentioned interpretation 
issue concerning the presence and origin of the BE distribution factor in the FDT. 
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4 Comments and conclusions 



First of all, we point out that, in order to construct the above mapping, the key ingredient 
we used is the hypothesis that linear response theory is applicable, which is the main 
hypothesis under which the FDT is established. When this is not the case, Eq.dlj) cannot 
be derived and we do not arrive to Eqs. (ITOl) and (ITTl) . which are crucial to build up our 
mapping. 

Note now that, by considering the "equivalent" harmonic oscillators system Sosc rather 
than the real one, we are somehow allowed to regard the BE distribution factor coth (^^) 

of the FDT in Eq. ([1]) as originating from the individual contributions coth ^ ^^j' j of each 

of the oscillators of the fictitious system (see above, Eqs. (129 p . (130|) and ( 132|) ). In this sense, 
such a mapping allows for an oscillator interpretation of the BE term in the FDT. 

At the same time, however, our result shows that this BE factor does not describe 
the physics of the system, i.e. it does not encode any real, physical, harmonic oscillator 
degrees of freedom of the system (see also the considerations below). 

In this respect, it is worth to point out that what we have implemented is not a 
canonical transformation, i.e. it is not a transformation which allows to describe the system 
in terms of new degrees of freedom (such as normal modes), but a formal mapping, a 
mathematical construct, which can be established, we repeat ourselves, only within the 
framework of linear response theory. 

In our opinion, then, our finding provides an answer to the questions of the "physical 
meaning" or "physical origin" of the BE term in the FDT or, stated differently, to the 
question of whether this BE distribution factor possibly describes the physical nature of 
the system or not . 

In fact, from the derivation of the FDT, we know that the BE factor derives from 
a peculiar combination of Boltzmann factors (see [7] and Eq. ([9]) above). At the same 
time, we have shown that, regardless the bosonic or fermionic nature of the (real) system 
S, it is always possible to establish a mapping which relates iS to a system of harmonic 
oscillators Sosc so that this BE factor can be regarded as "originating" from the individual 
oscillators of the "equivalent" system Sosc- Therefore, it is not the physical nature of the 
system which is encoded in this BE term but rather a fundamental quantum property of 
any bosonic and/or fermionic system: whenever linear response theory is applicable, any 
generic system is, at least with respect to the FDT, equivalent (in the sense defined above) 
to a system of quantum harmonic oscillators. 

Before ending this section, it is probably worth to spend few words on some examples 
of realistic systems where our mapping is at work. In this respect, we would like to 
note that there are several applications in the literature where fermionic systems, after 
bosonization, are actually described by a system in interaction with a bosonic bath. This 
is, for instance, the case of the anisotropic Kondo model, which is shown to be equivalent 
to a spin-boson model (a two level system in interaction with a bosonic bath). The same 
is also true for a quantum dot interacting with external leads. 

In the above examples, the system is described with the help of a spin-boson Hamil- 
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tonian, thus providing concrete realizations of the mapping discussed in this work. What 
amounts to the same thing, they are worked out examples where the Caldeira-Legget 
model is explicitly derived. 

In this respect, in fact, it is important to note that any application of our mapping is 
concretely substantiated in the Caldeira-Leggett model. At the same time, we stress again 
that the present work is focused on the question of deeply understanding what is really 
behind the fact that this modelization is so successful in covering the essential features 
of dissipative systems. We believe that we achieve this goal by performing a thorough 
analysis of the fluctuation-dissipation theorem. 

In summary, we have found that when linear response theory applies, any generic 
system can be mapped onto a fictitious system of harmonic oscillators so that that the 
mean square (A^(co')) of the fluctuating observable and the corresponding imaginary part 
of the generalized susceptibility Xa ('^) ^^^^ system are given by the corresponding 

quantities of the fictitious one. Moreover, we have seen that such a mapping allows to 
consider the BE distribution factor which appears in the FDT as originating from the 
individual harmonic oscillators of the fictitious equivalent system. This strongly suggests 
that it is only in this sense that this BE factor can be interpreted in terms of harmonic 
oscillators and that no other physical meaning can be superimposed to it. 

We believe that our mapping has a broader range of applicability than the worked 
case of the FDT discussed in this paper. Work is in progress in this direction. 
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